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The directed polymerization of dendritic actin networks is an essential element of many biological 
processes, including cell migration. Different theoretical models for the interplay between the un- 
derlying processes of polymerization, capping and branching have resulted in conflicting predictions. 
One of the main reasons for this discrepancy is the assumption of a branching reaction which is first 
(autocatalytic) versus zeroth order in filament density. Here we introduce a unifying framework from 
which these two models emerge as limiting cases for low and high filament density, respectively. A 
smooth transition between these cases is found at intermediate conditions. We also derive a thresh- 
old for the capping rate, below which zeroth order characteristics are predicted to dominate the 
dynamics of the network for all accessible filament densities. For capping rates above this threshold, 
autocatalytic growth is predicted at sufficiently low filament density. 



In many situations of high biological relevance, in- 
cluding the migration of animal cells and the movement 
of intracellular pathogens such as the bacterium Liste- 
ria, motility and force generation results from the di- 
rected polymerization of a dendritic actin filament net- 
work The organization of the growing network is 
determined mainly at the leading edge, where a small 
number of proteins regulates the interplay between three 
fundamental processes. The driving force for propulsion 
is polymerization of actin filaments from globular actin 
monomers, which is directly counterbalanced by filament 
capping through specialized capping proteins. In addi- 
tion, daughter filaments are branched off from mother 
filaments 0]. Although the biochemical details of this 
process are not completely understood, it is widely ac- 
cepted that the branching complex Arp2/3 is activated 
by nucleation promoting factors (NPFs). When an acti- 
vated Arp2/3-complex is bound to the side of an existing 
actin filament, a daughter filament starts to grow at a 
characteristic angle around 70° relative to the mother 
filament (compare Fig. [T^,). The interplay between the 
three fundamental processes of polymerization, capping 
and branching is force-dependent and leads to a non- 
trivial force-velocity relation that is essential for under- 
standing e.g. cell motility. 

Due to the high biological relevance and universal na- 
ture of the underlying processes, many theoretical models 
have been suggested to describe the characteristic fea- 
tures of growing actin networks However, in many 
cases contradictory predictions have been obtained, in 
particular regarding the force- velocity relation|3-yjJ and 
the spatial organization of the network [llHl3 |. Inter- 
estingly, many of these contradictions are a direct conse- 
quence of two different choices for the order of the branch- 
ing reaction. In autocatalytic models, the total rate of 



branching is assumed to be proportional to the num- 
ber of existing filaments in the network, i.e. it is mod- 
eled as a first order reaction in filament density, implic- 
itly assuming an unlimited reservoir of activated Arp2/3 
[H EEGll. Thi s yields growing actin networks for which 
a constant filament density is maintained only at a unique 
steady state growth velocity. Increasing forces acting 
against the network reduce the speed of growth only tran- 
siently, as an increasing filament density subsequently 
lowers the force per filament back to the stationary level. 

In marked contrast to the autocatalytic models, an- 
other class of models assumes that branching occurs with 
a constant rate, i.e. it is taken to be a zeroth order reac- 
tion in filament density, corresponding to a limited sup- 
ply of activated Arp2/3 [H, [TEES- Under these condi- 
tions, it has been shown that a continuum of steady state 
velocities exists. Moreover two competing steady state 
filament orientation patterns are stable, namely the ±35 
and +70/0/— 70 patterns shown schematically in Fig. [TJb 
and c, respectively. Transitions between these two fun- 
damentally different network architectures can be trig- 
gered by changes in network growth velocity [l7| . Indeed 
similar structural transitions have been demonstrated re- 
cently with electron microscopy for the lamellipodium of 
keratocytes [HI, EH. In this Letter, we will show that the 
two contradictory model scenarios of autocatalytic and 
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FIG. 1: (a) Interplay of polymerization, capping and branch- 
ing at the leading edge of an actin network growing towards 
the top. (b) A ±35 pattern is usually associated with den- 
dritic actin networks, (c) However, for certain parameter val- 
ues a +70/0/— 70 pattern can become stable. 
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zcroth order branching can be unified within a general 
theoretical framework that reconciles some of the seem- 
ingly contradictory observations and predictions. 

Model definition. Because the lamellipodium of cells 
migrating on flat substrates is very thin, we consider a 
two-dimensional model for a growing actin network. Fol- 
lowing earlier work along these lines [Tl[ [ID, [I?} , we con- 
sider the distribution function N(9, t) for the number of 
uncapped filaments per area orientated at time t at the 
angle 9 with regard to the normal of the leading edge. We 
assume N(9, t) to be defined in a small reaction region 
of width dbi- behind the leading edge, in which Arp2/3 
activation by NPFs and subsequent branching of new fil- 
aments is possible. Its time evolution is described by the 
kinetic equation 



dN(9,t) 
dt 



J W(0M')N(0' ,t)de' 
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(1) 



-k c N(e,t)-kUv nw )N(e,t), 



where the three terms on the right hand side describe the 
effects of branching, capping and outgrowth from the re- 
action region, respectively While capping is assumed to 
be a first order reaction in TV with rate k c , the order of 
the branching reaction is assumed to be variable and de- 
scribed by the dimensionless parameter \i. The two cases 
jjL = and /i = 1 correspond to zeroth and first order 
branching, respectively. Values of < \i < 1 interpolate 
between these two cases, fcb is the rate of branching and 
W(0, 9') = W{\0-8'\) is the distribution function for the 
orientation 9 of a daughter filament branching off from 
a mother filament with orientation 9' . Motivated by ex- 
perimental observations, we approximate this function by 
the renormalized sum of two Gaussians centered around 
±70° and each with standard deviation 5° (the detailed 
values are not relevant for our results). The third term 
in Eq. ([1]) describes the rate at which filaments are grow- 
ing out of the reaction region. We assume a well-defined 
network growth velocity v nw , which has to be compared 
with the single filament polymerization velocity vsx- For 
\9\ < arccos(f nw /wfii), the single filaments can keep up 
with the leading edge and we have fc| r (w nw ) = 0. If the 
orientation angle exceeds the threshold, however, the sin- 
gle filaments grow too slowly and leave the reaction re- 
gion with a rate fcg r (w nw ) = (v nw — vm cos 9) / <i br ■ 

Our model has only four relevant parameters, the rates 
k c and k^ for capping and branching, respectively, the 
network growth velocity v nw and the order fi of the 
branching reaction. We solve it numerically by discretiz- 
ing Eq. ([I} in 360 values for 9. In the following, we com- 
pare these results to two other versions of the model [13] ■ 
In a more detailed model, which does not use the distri- 
bution function N(9,t), we perform stochastic computer 
simulations of individual filaments evolving according to 
the rules described above. In a simplified analytical ver- 
sion of the reaction model Eq. (fTJ), we integrate over 35° 
sized angle bins and neglect contributions of filaments 



with growth directions > 87.5°. This leads to a set of 
three coupled equations: 



dAf O° _ 1 U. N ±70° _ jL A/no 



-(k c + k 3 J")N ± 

X - I &c + k lf) N ±70 



(2) 



dN ±m o _ ^ JV, 



dt 



where the summation is performed over the angle bins 
i = {0°,±35°,±70°}. 

Steady state analysis of the analytical model. From 
Eq. @, one can show that our model has two steady 
state solutions, 7V ss7 ° and iV ss35 : 
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and 



where 



A = k h /^2k c (k c + k™°), B = fc b /2(fc c + k 3 f). (5) 

The two fixed points of the system correspond to the two 
competing orientation patterns depicted schematically in 
Fig. [TJi and b, respectively. While iV ss70 has peaks at 
+70/0/-70 degrees, 7V ss35 is peaked at ±35 degrees. 

Linear stability analysis shows that for (j, > 1, the two 
fixed points are both saddle points and thus no stable so- 
lutions exist. For fx < 1, the eigenvalues for the two fixed 
points have opposite signs, which leads to mutual exclu- 
sive stability. Fig. [5^ illustrates the regions of stability 
for each of the two orientation patterns within the three 
dimensional parameter space spanned by fct>, &c and u nw . 
The two meshed surfaces represent transitions between a 
+70/0/— 70 pattern above and below the two surfaces and 
a ±35 pattern in between. These transition surfaces are 
independent of [i and fcb ■ Thus all cases with fj, < 1 show 
no difference compared to a model explicitly assuming 
zeroth order branching [l7j |. 

The case = 1 (autocatalytic growth) is special and 
has to be discussed separately. In the limit /j, — > 1, each of 
the two steady state solutions can exist only if additional 
conditions on the parameter values are satisfied. This 
important aspect can be understood from Eq. |T]), where 
for /i = 1 all three terms on the right hand side become 
linear in N and thus for the steady state can sum up 
to zero only if the three rates cancel each other. The 
corresponding condition is different for each of the two 
fixed points and can be inferred from Eq. ([3|) and Eq. (j4} , 
respectively, in combination with Eq. ([S]). For N ssJ0 , only 
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FIG. 2: Phase diagram of the analytical model. In general, two mutually exclusive stable fixed points exist, the ±35 and 
+70/0/— 70 patterns shown in Fig. la and b. (a) For /i < 1 the parameter space is three-dimensional, with the ±35 pattern 
being stable between the two meshed surfaces. The ±70/0/— 70 pattern is stable above and below these two surfaces. For the 
case fj, = 1 (autocatalytic growth), only one network growth velocity is possible and the stable solutions are restricted to the 
dark and light gray surfaces for the ±35 and ±70/0/— 70 patterns, respectively. The simplifications of the analytic model lead to 
a jump discontinuity at the lower transition, (b) Autocatalytic model (/i = 1) in the fc c -fcb plane. Here the jump discontinuity 
is not visible. The inset demonstrates the good agreement between the analytical and the full numerical solutions. fc™ m and 
£;™ ax show the parameter region over which autocatalytic growth is possible for the given effective branching rate k^ c (see 
Eq. (|10p and Eq. below), (c) Autocatalytic model in the k c -v n w plane. Here the jump discontinuity leads to a region 
where no velocity with finite N exists. The result from the full numerical solution (inset) does not share this property. 



A = 1 gives a finite solution for fj, — > 1, thus we require 
klf{v nw )k c ± k 2 c = kl/2. For N ss35 , only B = 1 gives 
a finite solution for \i — > 1, thus we require (v nw ) ± 
k c = kb/2. We conclude that for fi = 1, a finite filament 
density in steady state is only observed at one unique 
network velocity for each combination of kb and k c , as is 
expected for autocatalytic growth fl5| . 

In Fig. the two conditions resulting from A = 1 and 
B = 1 arc shown as light and dark gray surfaces, respec- 
tively. Because the simplified model is too restricting as 
to be able to realize any value for network growth veloc- 
ity, there is a jump discontinuity to the lower transition 
(at the upper transition, the two surfaces coincide). In 
Fig. [2Jd and c, we show the projections of the phase dia- 
gram for (i = 1 from Fig. [5^ onto the k c —kb and k c -v nw 
planes, respectively. In Fig. Oj, the jump discontinuity 
is not visible. In Fig. it leads to a region with ve- 
locity values that do not yield finite Ns\. Note however 
that this special property only occurs as an artefact of 
the simplications made for analytical progress. The in- 
sets in Fig. [5b and c show the phase diagrams obtained 
from the full numerical analysis of Eq. (JTJ) . In contrast 
to the analytical model, they show no jump discontinuity 
and a larger stability region of the ±35 pattern at high 
capping rates. Otherwise the agreement is very good, 
demonstrating the strength of the analytical results. 

Reaction model for branching. Until now we have 
treated the order [i of the branching reaction as a phe- 
nomenological constant. In order to relate (i to filament 
density, we now introduce a reaction model for branch- 
ing, based on a likely scenario for the involved biochem- 
ical steps [HI, El- We consider two variables: Arp is the 
concentration of Arp2/3 that is bound to the filaments, 



but did not lead to a daughter branch yet. Npf is a frac- 
tion of the total concentration of obstacle bound NPFs 
at the leading edge (Npfo) which are available to activate 
bound Arp2/3 complexes to nucleate a daughter branch. 
The kinetic equations are 



dArp 

dt 
dNpf 

dt 



k+N m 



Arp - fc b Arp Npf , 



-fcbArpNpf ±fc act (Npf -Npf) 



(6) 



Arp increases as more complexes bind to the filaments 
with rate fc+. It decreases due to dissociation (rate fe_), 
outgrowth (rate iw/^br) and activation (rate kb)- This 
activation step also decreases Npf as we assume NPFs 
that activate Arp2/3 to be occupied for additional inter- 
actions with other Arp2/3 complexes until they become 
available again at rate fc ac t- 

In steady state, the effective rate of branching is 



BR SS (JV m ) = fc b Ar Pss Npf ss . 



(7) 



which is shown in Fig. [3] for a typical set of parameters. 
At low filament number, we can expand the right hand 
side of Eq. ([7]) to linear order in Nfu: 



BR S Q S 



Npfodbr^+^b 



Npfo dbr kb + dbrk- ± v 



-JVfii = K c N m (8) 



thus defining an autocatalytic rate k^ . In the limit of 
large filament density, BR SS approaches a constant: 



BR S ^ = Npf fc a 



(9) 



and hence a zeroth order branching rate. Thus the cases 
of autocatalytic growth (/i = 1) and zeroth order branch- 
ing (/i = 0) emerge as the limits for low and high filament 
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FIG. 3: Effective branching rate BR SS as a function of the 
number of filaments Na\ in steady state growth. For low fila- 
ment number the branching reaction is linear (autocatalytic) 
as given in Eq. (0 (dashed line), while in the limit of high 
filament number zeroth order branching is observed with a 
constant rate given by Eq. (dotted line). The inset shows 
the corresponding steady state concentrations of Arp (solid 
line) and Npf (dashed line). 



density, respectively. At intermediate filament densities 
a smooth transition between these two limiting cases oc- 
curs. 

We next couple the two models Eq. (p} and Eq. (|6|). For 
steady state growth, this can be done by either deduc- 
ing the reaction order parameter fi as a function of Nfn 
from a comparison of the branching rates in both mod- 
els, or by directly inserting a filament density dependent 
branching rate, fcb = BR^A^i), in the branching term 
of an otherwise zeroth order model with /.t = 0. The 
second strategy turns out to be numerically more robust 
and will be used in the following for calculating station- 
ary filament densities at given network growth velocity 
in steady state. 

Limits of autocatalytic network growth. For a growing 
network at low filament density, well within the first order 
(autocatalytic) regime, an effective first order branching 
rate can be approximated by the coefficient of the linear 
expansion, k^ c in Eq. (|8]). Then [i = 1 and the contribu- 
tions from capping and outgrowth of filaments integrated 
over all filament directions have to exactly balance fila- 
ment nuclcation via branching, compare Eq. ([1]). From 
this condition it is possible to derive fundamental lim- 
its of autocatalytic growth in terms of minimum and 
maximum capping rates k™ m and fc™ ax , corresponding 
to the largest and smallest possible outgrowth rate at 
iw/^fii = 1 and v nw /vfn = 0, respectively. By inserting 
the autocatalytic branching rate fc ac defined in Eq. © 



FIG. 4: Network growth velocity as a function of filament 
density for different capping rates obtained from the numeri- 
cal solution of the full model with filament density dependent 
branching rates (full lines). For comparison, we show the ze- 
roth order description with constant branching rate (dashed 
lines). For capping rates fc™ m < k c < fc™ ax , an autocatalytic 
regime with is observed at low filament density. Here the 
solid plateau does not agree with the dashed lines. A capping 
rate k c = fc™ ln (thick black solid line) marks the transition 
to a zeroth order behavior, when the solid and dashed lines 
start to agree. The inset shows the results from a stochastic 
computer simulation, which are in excellent agreement with 
the results from the rate equation model for the orientation 
distribution function. 



into the first order condition, we get 



^[(cos70°-l) 



and 



,(3!^- + ( cos70 o_;l)2 



^max b 



V2 



(10) 



[11) 



These two threshold values arc also illustrated in Fig. [5}}. 

Fig. [4] shows the steady states of a network at various 
capping rate k c by numerically solving the full model with 
an effective branching rate as a function of filament den- 
sity (solid lines) . For a capping rate k c in between the two 
limits Eq. (fTO)) and Eq. (fTT|) . an autocatalytic regime is 
observed at low filament density. Here, the network speed 
shows a plateau and thus several values for the filament 
density correspond to the same velocity, as it is typical 
for an autocatalytically growing actin network. Increas- 
ing network density above the intermediate regime, the 
steady state network velocity decreases like in a purely 
zeroth order description (dashed lines). The steady state 
network velocity in the autocatalytic domain drops as k c 
increases, until at around k c ~ fc™ ax the filament density 
exponentially decays to zero for all accessible network ve- 
locities. For a capping rate k c < fc™ ln , even the fastest 
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network velocity, tw/^fii = 1, is not able to balance fila- 
ment branching and capping via outgrowth anymore. In 
a purely autocatalytic description, this yields diverging 
filament densities. On the contrary, within our unifying 
framework the filament density is only able to increase 
until the transition to zeroth order growth takes place 
and thus a reasonable steady state is obtained at finite 
filament density. As a direct consequence, this yields Nri 
versus v nw curves (solid) with no significant difference 
to a purely zeroth order branching model (dashed) for 
such low capping rates. The inset shows the results of 
the stochastic network simulations, which are in excel- 
lent agreement with the numerical solution of the rate 
equation. This suggests that our results are very robust 
in regard to the details of the model implementation. 

Conclusion. In this Letter, we have shown that con- 
flicting predictions of two important classes of actin net- 
work growth models, namely first order branching (au- 
tocatalytic) and zeroth order branching models, can be 
resolved within a single theoretical framework. This uni- 
fying model yields first and zeroth order characteristics in 



the limits of low and high filament density, respectively. 
Using this framework we have also shown that a lower 
threshold for the capping rate exists, below which zeroth 
order characteristics are predicted to dominate network 
growth, thus avoiding unrealistic divergences in the fila- 
ment density. Although the exact activation cascade of 
Arp2/3 is a subject of current research, our arguments 
apply quite generally as long as network growth is lim- 
ited by the availability of activated Arp2/3. Our results 
suggest to conduct experiments which systematically ex- 
plore the different predicted growth regimes. 
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